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^^ ' Abstract 

o: 

rvj ^ The major role type la supernovae play in many fields of astrophysics and 

, . in particular in cosmological distance determinations calls for self-consistent 

O ' models of these events. Since their mechanism is believed to crucially depend 

.^^ , on phenomena that are inherently three-dimensional, self-consistent numerical 

• models of type la supernovae must be multi-dimensional. This field has recently 

>0 ' seen a rapid development, which is reviewed in this article. The different mod- 

^~~^ , eling approaches are discussed and as an illustration a particular explosion 

^_^ ' model - the deflagration model - in a specific numerical implementation is pre- 

ri \ sented in greater detail. On this exemplary case, the procedure of validating the 

Oj. model on the basis of comparison with observations is discussed as well as its 

' ' application to study questions arising from type la supernova cosmology. 

^ . 1 Introduction 

The fact that many astrophysical processes are inherently three-dimensional makes 

^ . realistic numerical simulations challenging. Although computational resources in- 

\Q ' crease steadily, complex phenomena cannot be directly resolved in such simulations 

VO > in the foreseeable future. Even exploiting the available computational power and 

'^ ' memory, as most multi-dimensional astrophysical simulations do, considerable ef- 

^ ; fort is required in modeling the processes in a way that despite the shortcomings 

■^ in resolution the results provide credible physical approximations to the problems. 

^^ . Moreover, the implementation of efficient numerical techniques is one of the corner- 
stones of successful multi-dimensional modeling of astrophysical phenomena. 

Type la supernovae (SNe la henceforth) are an excellent example for this class 
of astrophysical objects. Numerical simulations of these events in three dimensions 

\^ easily reach the limits of today's computational resources. But as the underlying 



o 
> 



?H physical mechanism is believed to crucially depend on three-dimensional phenom- 
. . ena, such as turbulence, this is the only con ceivable way towards self - consis tent mod- 



els of type la supernova explosions. Since lHillebrandt & Niemeyeg (120001) reviewed 



type la supernova explosion models, the field has witnessed a brisk development. 
Progress in numerical methods and computational capabilities facilitated pioneering 
three-dimensional (3D) simulations, which will be reviewed in the following. 



Self-consistent SN la models are called for in order to achieve a sound under- 
standing of the mechanism of these astrophysical events. This is motivated by their 
significant impact on many aspects of astrophysics and cosmology. Being one of the 
main sources of iron group elements, SNe la contribute to the chemical evolution 
of galaxies (e.g. Francois et al., 2004) . They affect star formation and drive shock 
waves in the interstellar and intergalactic media. 

Most remarkable, however, was the application of SNe la in observational cos- 
mology, where these object s were employed as distance indicators (as put forward by 



Branch & Tammannlll992h . Evidently, SNe la are a valuable tool to e xtend the Hub- 



ble diagram to large redshifts and to determine the Hubble constant jHamuv et al 



Il996; Branch. 1998) . A t redshifts above 0.5, a significant deviati on from the linear 
Hubble law was noticed (iRiess et al.[ll998HPerlmutter et all[l999l) . Here, SNe la ap- 
pear dimmer than expected in a matter-dominated flat or open Friedmann-Robertson- 
Walker Universe. This led to the spectacular interpretation that the expansion of 
Universe is currently undergoing an acceleration. The determination of the driving 
force of this acceleration is one of the greatest challenges in contemporary physics. 
Meanwhile it is parametrized as "dark energy" (see e.g. Leibundgut, 2001 ). Determi- 
nations of cosmologic al parameters based o n anisotropies in the cosmic microwave 
background radiation f Sperge l et al.L l2003b and on large-scale galaxy surveys pro- 
vided independent confirmation of the SN la measurements. 

Yet the question of the applicability of SNe la as distance indicators is still not 
satisfactorily answered. SNe la are remarkably uniform events by astrophysical stan- 
dards, but evidently no standard candles. Only a calibration of the distance measure- 
ments according to empirical correlations between observables provides the neces- 
sary accuracy for the determination of cosmological parameters. A firm theoretical 
reasoning of such correlations is, however, still lacking. 

The simplest form of dark energy is a cosmological constant, but more com- 
plicated contributions to the energy-momentum tensor in the Einstein equations are 
also conceivable. How can one determine the nature of the dark energy? A first step 
would be to constrain its equation of stat e. SNe la seem to be a suitable tool for this 
task and currently two major campaigns JAstier et al.l 120051 : ISoUerman et al.l |2005 j 
apply them in distance determinations of hundreds of supernovae out to redshifts 
of z ~ 1 and systematic satellite-borne SN la observations are planned. The large 
number of observations is necessary to reduce the statistical errors because putting 
tight constraints on dark energy equation of state (see Astier et al.. 2005, for recent 
results) requires a high accuracy of the distance determinations. But still, potential 
systematical errors arising from calibration of the distance measurements may ob- 
scure the results. Getting a handle on these is one of the goals of modeling SNe la. 
Obviously, the predictive power of self-consistent multi-dimensional SN la models 
is promising for progress in answering the questions arising from SN la cosmology. 

In Sect.|2]we set out how the astrophysical scenario of SNe la is derived from ob- 
servations. Numerical approaches to modeling SNe la are discussed in Sect. [3]with 
special emphasis on multi-dimensional models. As an illustrative example, a partic- 
ular implementation of the deflagration SN la model is presented in Sect. H] where 
also ways to test the validity of SN la models on the basis of comparison with obser- 
vations are outlined and the application to questions arising from the cosmological 



application of SNe la is discussed. Conclusions are drawn in Sect.|5j 



2 Observations and astrophysical scenario 

Observational features of SNe la suggest a specific astrophysical scenario. The cor- 
nerstones of the astrophysical model of SNe la are set by two fundamental char- 
acteristics of these events. Evidently, SNe la belong to the most energetic cosmic 
explosions, releasing about 10^^ erg of energy. For a short period of time they can 
outshine an entire galaxy consisting of tens of billions of stars. SNe la spectra are 
characterized by the lack of indications for hydrogen and helium which together 
with a pronounced P Cygni silicon line at maximum light classifies these objects 



dWheeler & Harknesslll99 0!). Lines of intermediate mass elements (such as Si, Ca, 



Mg and S) and oxygen are observed in near-maximum light spectra (e.g. lFilippenkol 



1997blal) . With respect to light curve and spectra observatio ns SNe la form a class 



of remarkable homogeneity (e.g. Branch & TammannLll992 l) 



Assuming supernovae to originate from single stellar objects, only th eir gravi 



tation al binding energy, released in a collapse towards a c ompact object jZwicky , 



tational bmamg energy, releasea m a collapse towaras a c ompact oDiect QZ-wickv , 
19381) . or its nuclear energy, released in explosive reactions fHo yle & Fowlen.ll960r) . 



come into consideration as possible energy sources. In the particular case of SNe la 
no compact object is found in the remnant excluding the first possibility. The homo- 
geneity of the class of SNe la and the fact that no hydrogen is found in their spectra 
provides a strong hint that the object undergoing the nuclear explosion may be a 
white dwarf (WD) star consisting of carbon and oxygen (Ch-O). 

Lightcurves of SNe la rise over a time scale of several days and decline over 
months. It is therefore clear that they cannot be powered directly by the explosion 
since the temperatures fall off mu ch too rapidly in the expansion. This problem was 
solved by' Truran et all 11967') and ' Colgate & McKeeJ (|l969i) who suggested that the 



^^Ni produced in large amounts in the explosive thermonuclear burning provides the 
energy source for the optical event by radioactive decay to ''^Co and ''^Fe. 

2.1 Progenitor evolution and ignition 

A single WD is an inert object. How can it reach an explosive state? The only way 
to introduce the necessary dynamics into the system is to assume it to be part of a 
binary system and to gain matter from the companion. Several models have been 
proposed for this progenitor evolution. 



In the double degenerate scenario (llben & Tutukovill984tlWebbinklll984i) . two 



C-hO WDs merge. The lighter of the two is disrupted and its matter accreted onto the 
heavier WD. To become a potential candidate for a SN la, the merger should even- 
tually reach the Chandrasekhar-mass (i.e. of the maximum mass that is supported 
against gravitational collapse by the pressure of the degenerate electrons). Stellar 
evolution predicts that systems fulfilling this requirement should exist and indeed a 
potential candidate that w ill merge in less than one Hubble time has been detected 



(iNapiwotzki et al.L l2005h . Mergers of two C+0 WDs provide a natural explanation 



for the absence of hydrogen. However, numerical simulations indicate that the high 



accretion rate onto the more massive WD leads to an off-center ignition and the sub- 
sequent burning could convert the material to oxygen, neon, and magnesium. A WD 
of this composition, however, tends to undergo a gravitation al collapse rather than a 
thermonuclear explosion (e.g. lSaio &Nomotolll985ill998h . 



In the sing l e deg enerate scenario (IWhelan & Ibeni Il973t iNomotol Il982 : 



Iben & Tutukovi 119841) . the WD accretes matter from a non-degenerate compan- 
ion (either a main sequence or an AGB star). This idea was recently supported 
by the detection of t he po tential companion of Tycho Brahe's 1572 supernova 



(IRuiz-Lapuente et al.L 120041) . which is a solar-type star. A detailed analysis of this 



scenario, however, revealed that the accretion rates admissible here are restricted to 
a narrow window. Too low rates evoke to nova eruptions in which the WD looses 
more material than accreted before and too high rates would lead to the formation 
of an extended He-rich envelope. Moderate accretion rates build up a degenerate 
He-shell which could detonate and trigger a detonation of the carbon-oxygen core. 
Since this happens before the WD reaches the Chand rasekhar-mass limit, such a 
SN la model is termed sub-Chandmsekhar explosion ( Wooslev & Weaveg, 119941) . 
Somewhat higher accretion rates, however, can lead to quiet hydrostatic burning of 
the accreted material processing it to carbon and oxygen. In this case, the WD may 
reach the Chandrasekhar mass {Chandra sekhar-mass mot/e/. iHoyle & Fowlen.ll960t 
lArnett. 1969 : Hansen & Whee lej. 19691). Potential candidates fo r such systems are 
Supersoft X-ray Sources (e.g. lKahabka & van den HeuvelLll997l) The limitation the 
fuel available in the explosion to the Chandrasekhar mass (~1.4M0) makes the de- 
scribed scenario particularly favorable since it provides a natural explanation for the 
striking uniformity of SNe la in the gross observational features. On the other hand, 
it is afflicted with great uncertainties. Achieving a stable mass transfer in the progen- 
itor binary system to b uild up a Chandrasekhar mass WD is highly non-trivial (e.g. 
Nomoto & Ibenlll985h and the observational evidence for such systems is sparse. 



Although the striking homogeneity holds for most SNe la (so-called Branch- 
normals, £rancheLaLJ^93), s ome events differ significantly. These explode rn uch 
weaker (like SN 1991bg, see iFihppenko et alj [l 992a; Leibun dgut et al J 1 1 993h or 
more vigorous (hke SN 1 991T, see |P hillips et al. 1992; Filippen koet alJll992bt 
Ruiz-Lapuente et al.lll992t IJeffervetaLi 1992 : Spvromilio et al. 19921) than the av- 
erage. At present it is unclear, whether the entire class of SNe la can be explained by 
only one progenitor scenario. The currently favored Chandrasekhar-mass model may 
possibly only be able to explain the Branch-normals while sub- and superluminous 
events require different progenitor scenarios. Therefore, even if currently not in the 
focus of research, the double-degenerate and the sub-Chandrasekhar mass scenarios 
may contribute to the SN la population. In the following, however, we will focus on 
the Chandrasekhar-mass model, since it received most attention in recent theoretical 
modeling. 



When the WD approaches the Chandrasekhar limit, the density at the center of 
the WD increases rapidly so that fusion of carbon ignites. Contrary to the situation 
in main sequence stars, the degenerate material of the WD does not allow for mod- 
eration of the burning by expansion. Heat transport is achieved here by convection 
giving rise to a stage of convective carbon burning that lasts for several hundred 
years. This phase is terminated by one or more small spatial regions undergoing a 
thermonuclear runaway, marking the birth of a thermonuclear flame and the onset 
of the explosion. The convective burning stage and the conditions at flame ignition 
are extremely hard to model both analytically and numerically. Therefore the ex- 
act shape and location of the first flame spark(s) is not yet well constrained. These, 
however, are crucial initial parameters in multi-dimensiona l explosion models. Only 



a few studies addressed the flame ignition process so far iGarcia-Senz & Wooslev 



( 119951) simulated the thermonuclear runaway of a hot bubble floating upward from 



the center of the WD. They concluded that off-cente r and multi-spot ignitio ns are 



possible. A simular result was recently obtained by llapichino et alJ (120051) . The 
convective phase directly preceding the ignition crucially influences the configura- 
tion of the initial flame. Here different studies led to controversial results. While 



Wooslev et al.l (120041) and iKuhlen et al.l (12005[) favor off-center, po ssibly multiple 



and asymmetrically distributed flame sparks, iHoflich & SteinI (120021) put forward a 
single central ignition. 



2.2 Flame propagation and explosion 

The goal of SN la explosion models is to follow the propagation of the thermonuclear 
flame from its ignition near the center of the WD outwards and to determine the 
composition and the distribution of the burning products in the ejected material. 

Theoretically, SN la explosions are governed by the equations of reactive fluid 
dynamics, i.e. the Navier-Stokes equations extended with an equation of species bal- 
ance, suitable terms for heat conduction, diffusion of species, a nd source terms of en - 



ergy and species in combination with an equation of state (eg. lOran & Borislll987h . 
These equations allow for solutions of traveling reaction waves converting unburnt 
to burnt material. If the scales under consideration are much larger than the internal 
structures of these waves, a simplified description of the system is possible. Neglect- 
ing all phenomena that actually govern the propagation of the reaction wave, i.e. the 
reaction kinetics and transport processes, the burning front is modeled as a moving 
discontinuity. This picture provides a description of the hydrodynamical state of the 
material in the unburnt and burnt regions and the equations simplify to the reactive 
Euler equations. From these, the Rankine-Hugoniot jump conditions for the state 
variables accross the burning front follow (see Landau & Lifshitz. 1959) which al- 
low for two different modes of front propagation. One is the subsonic deflagration in 
which the flame is mediated by the thermal conduction of the degenerate electron gas 
and the other is a supersonic detonation in which the burning front is driven by shock 
waves. Either one of these modes or a combination of both have been suggested in 
different explosion models. 



The p rompt detonation model was first applied in a numerical simulation by 



Arnett (1969) . A spherically symmetric detonation wave was initiated near the 



center of the WD and propagated outward. Criteria for the spontaneous for 



mation of a detonation wave have been investigated bv lBlinnikov & Khokhlov 
( 11986) and Woosley ( 1990) . They reach the conclusion, that a detonation can 
initiate only under certain prerequisites. A prompt detonation as explosion 
model produces enough energy for a SN la event. However, ahead of a su- 
personic detonation wave the fuel cannot expand and is therefore incinerated 
at the high densities of an equilibrium white dwarf. This results in the al- 
most complete conversion of the material to nickel-peaked nuclear statistical 
equilibrium (Arnett, 1969 ; Arnett et al., 1971) , which is in conflict with the 
intermediate mass elements observed in SN la spectra. These nucleosynthetic 
problems rule out a pure detonation scenario as a standard model for SN la 
explosions. 



The deflagration model ( INomoto et al.L Il976h assumes the flame propagating 



in the subsonic deflagration mode. The laminar burning speed of the deflagra- 
tion flame is determined by microphysical transport proc esses. For conditions 



of ca rbon burning in C-i-O WDs it is highly subsonic dTimmes & Wooslev , 



1992h and therefore the flame propagates far too slowly to explain SN la ex- 
plosions. The expansion of the star will quench burning before the WD gets 
unbound. On the other hand, this model can cure the problem of nucleosyn- 
thesis, since rarefaction waves travel ahead of the flame with sound speed and 
lower the fuel density prior to burning. Thus the material can partly be pro- 
cessed into intermediate mass elements. 

The deflagration model undergoes a significant improvement when multidi- 
mensional effects are taken into account. The propagation of the deflagra - 
tion front is subject to several instabilities (e.g. Nie meyer & W oosley. 1997). 
Of purely hydrodynam ical origin is the Landau-Darrieus instability ( Landaui 



1944t iDarrieusL 119381) wh ich in the nonlinear stage is stabiUzed in a cellular 
pattern dZeFdovichi 119661) thus enlarging the flame surface area and enhanc- 



ing the net burning rate. The major effect accelerating the flame, however, 
is due to the buoyancy unstable flame propagation from the center of the star 
outwards. It leaves behind light and hot ashes below the dense fuel - a density 
statification inverse to the gravitational acceleration. In its non-linear stage, 
the Rayleigh-Taylor instability leads to the formation of mushroom-shaped 
burning bubbles raising into the fuel. The Reynolds number typical for this 
situation is as high as lO''*. Clearly, shear (Kelvin-Helmholtz) instabilities at 
the interfaces of these bubbles will generate turbulent eddies which then decay 
to smaller scales forming a turbulent energy cascade. The flame will interact 
with these eddies down to the Gibson-scale at which the turbulent velocity 
fluctuations become comparable to the laminar flame speed. Below the Gib- 
son scale, the flame burns faster through turbulent eddies than they can deform 
it, and the flame propagation is thus unaffected by turbulence there. This in- 
teraction corrugates the flame again increasing its surface and consequently 
accelerating the effective propagation speed. 



The delayed detonation (DD) model conjoins the advantages of the deflagra- 
tion and the deton ation models. It was put forwar d by Ivanova et al. (1974) . 



Khokhlovl ( Il991ah . and IWooslev & Weaver! (119941) . Burning starts out in the 



slow deflagration mode pre-expanding the star. Motivated by transitions from 
the deflagration to the supersonic detonation mode (deflagration-to-detonation 
transition, DDT) observed in terrestrial turbulent combustion processes, such 
a phenomenon is assumed to occur in SNe la. Its physical mechanism, how- 
ever, remains unclear and therefore it enters the model as a free parameter. 
Usually, in one-dimensional (ID) simulations the DDT is artificially initiated 
once the flame reaches fuel of a certain transition density pu-. The assumed 
detonation then burns the star until the flame is quenched by the expansion. 
This detonation is an easy way to explain the energy release necessary for a 
SN la explosion. The important notion in this model is that a detonation in low 
density fuel (pre-expanded in the deflagration stage) can lead to only partial 
burning and is therefore capable of generating intermediate mass elements. 
Another possible advantage of this model derives from a problem of current 
3D implementations of the pure deflagration scenario. The Rayleigh-Taylor 
bubbles being the origin of the turbulent flame acceleration cause consider- 
able amounts of unburnt C+O matter to remain in "fingers" near the center 
of the star, which are possibly in conflict with spectral SN la observations. 
A detonation wave initiated at later stages of the evolution could be capable 
to burning out those fingers and to process the previously unburnt material. 
The best agreement with observations was achiev ed for pa ~ 1 x lO^gcm"^ 
(IHoflich & Khokhlovl 1 1996tllwamoto et al.lll999h . The main disadvantage of 
the DD model is that this transition density remains an arbitrary parameter un- 
less the mec hanism of a possible DDT is physically determined (if there exists 



one at 



all,cf.|N 



lemever , 



1999h . 



The pulsational delayed detonation (FDD) model (^ Arnett & Livnel Il994allbh 
is similar to the DD model in the sense that it combines an initial deflagration 
with a later detonation. The flame is assumed to propagate in the initial defla- 
gration phase with its laminar burning speed and pre-expands the star. Due to 
the slow flame velocity, the burning front stalls and fails to unbind the star. The 
WD then re-contracts giving the interface between burnt and unburnt material 
enough time to mix and to become nearly isothermal. Compressional heating 
finally triggers a detonation at densities that are lower than that prior to the first 
expansion phase. iHoflich & Khokhlovl (119961) employ this scenario in a phe- 
nomenological ID model and conclude that it may account for sub-luminous 
SN la. However, the assumption that the flame propagates with the pure lam- 
inar burning velocity in the deflagration phase seems unrealistic, because of 
the flame instabilities and the resulting turbule nt flame acceleration. Recent 
multi dimensional deflagration models (Reineck eet al.[ l2002bt iGamezo et al. 



12003 ) demonstrated that taking into account these effects, the star is likely to 
get unbound instead of recontracting. 



3 Numerical models 

3.1 Relevant scales 

The numerical implementation of SN la models accounting for the full exploding 
WD star in multiple dimensions is significantly complicated by the wide range of 
relevant length scales involved in the problem. From the radius of the WD star 
(~2000 km at the onset of the explosion and expanding in the process) it reaches 
down to the flame width which is well below one centimeter. In the deflagration 
model the relevance of turbulent effects amplifies the scale problem since the turbu- 
lent cascade extends to the even much smaller Kolmogorov scale where the turbulent 
energy is dissipated into heat. Here, the flame interaction with the turbulent cascade 
down to the Gibson scale must additionally be taken into account. Current 3D sim- 
ulations capturing the entire star reach resolutions around one kilometer while the 
Gibson scale is of the order of 10"* cm at the beginning of the explosion and decreases 
steadily. 

For large-scale multi-dimensional SN la simulations this has three consequences. 

1 . The internal flame structure cannot be resolved. Thus, an effective flame model 
has to be applied and complementary small scale simulations are required. 

2. It is not possible to fully resolve the interaction of the flame with turbulence. 
Therefore modeling of the effects on unresolved scales is necessary. 

3. Assumptions about the flame properties at unresolved scales (e.g. stability be- 
low the Gibson scale) have to be validated in separate small-scale simulations. 

3.2 Modeling approaches 

Numerical models of SN la explosions have to face three major challenges. Apart 
from the vast range of relevant length scales they need to take into account inherently 
three-dimensional physical phenomena and to solve the kinetics of nuclear burning. 
To meet all these requirements in a single simulation will be impossible in the fore- 
seeable future. Therefore the problem has been tackled in different approaches. 

The first path towards SN la explosion modeling is to restrict the simulations to 
only one spatial dimension. Here, in principle a resolution of the relevant scales is 
achievable and a detailed description of the nuclear reactions is feasible. However, 
crucial three-dimensional physical mechanisms are not explicitly taken into account 
and need to be parametrized. Although such models were shown to yield good fits to 
observations and pioneered the insight into fundamentals of the explosion process, 
they lack consistency and are of little predictive power. 

In multi-dimensional simulations, contrariwise, the computational costs of mod- 
eling the explosion hydrodynamics is prohibitive to directly resolve all relevant 
scales as well as details of the nuclear processes. While the latter may be improved 
in the forthcoming years, and is meanwhile separated from the actual explosion sim- 
ulations still maintaining a reasonable accuracy (see below), even a drastic increase 
in computational capabilities will not allow for a resolution of all scales in multi- 
dimensional simulations. 



A third approach is to study specific effects on a limited range of spatial scales, 
to validate assumptions and improve modeling techniques of the large-scale SN la 
simulations. 



3.3 One-dimensional simulations 

Although we focus on multi-dimensional models of SN la explosions, we will first 
give a brief overview of ID spherically symmetric simulations. Since (at least in 
the initial deflagration stage) the propagation velocity of the thermonuclear flame is 
determined by multi-dimensional effects, such as instabilities and interaction with 
turbulent velocity fluctuations, a physically undetermined speed is ascribed to the 
spherical burning front. This free parameter can be used to fit the observations and 
therefore give a hint to the expected average temporal evolution of the flame in multi- 
dimensional s i mulat ions. The best known of the se models, W7, was presented by 
Nomoto et alj (119841) and Thielemann et al. (119861) and has been emp loyed in several 



studies since then (e.g. Iwamo to et al.Lll999tlBrachwitz et aUl20 00V Here the flame 
propagation speed was modeled with mixing length theory and started slowly to pre- 
expand the star. Later it strongly accelerated and consumed large fractions of the star 
converting the material to iron-group elements (mainly ^^Ni) and intermediate mass 
elements. Different ways to pa r ameterize the flame propagation veloci t y were pro- 
posed by iNomoto et all d 19761) . IWooslevI (Il990l) . iHoflich & Khokhlovl Jl996l) . and 
iNiemever & Wooslevl(ll997h . All these models conclude that a flame acceleration to 
one third of the sound speed is necessary for consistency of the results with observa- 
tions. The problem of overproduction of neutron-rich iron-group material noted in 
these models may be significantly reduced when taking into account revised electron 
capture rates (B rachwitz et al., 2000) . 

When applying ID models to simulate delayed-detonation scenarios, the con- 
dition for the DDT comes into play as a second undetermined parameter Such 
models are capable of providing excellent fits to lightcurves and s pectra and seem 
to yield reasonable compositions of the nucleosynthesis products ( Wooslevl Il990t 






Hoflich & Khokhlovl [l996l:lKhokhlovill991bl;llwamoto et al.i[l99 

Although lacking consistency and therefore predictive power, ID mod 
els have been employed to address questions arising from SN la cosmology , 
such as th e origin of the diversity in these objects (Bravo et al., 1993, 199' 
Hoflich & Khokhlov. 1996; Hoflich et al., 1998; Umeda et al., 1999; Iwamoto et al 



1999HDominguez & Hoflic"hll2000HDominguez et al.Ll2000ll200lh . 



3.4 Multi-dimensional simulations 

Initiated by iMiiller & Arn ett (^1982, 1986), multi-dimensional models were applied 
to fix the turbulent flame propagation velocity undetermined in spherically symmet- 
ric simulations. Given the wide range of scales on which the flame is affected by 
turbulence, this is an ambitious project, additionally challenged by the lack of reso- 
lution of the thermonuclear flame structure. For both problems, different approaches 
have been taken. All these were guided by the theory of turbulent combustion in 



terrestrial flames (see|Peters|,|200] 



3.4.1 Hydrodynamics 



The hydrodynamical equations are discretized either in a Lagrangean or Eule- 
rian approach. The Sm ooth Particle Hydrodynamics (SPH) technique used by 
Garcia-Senz et alj (Il998h falls i nto the first c ategory, while an implicit scheme on 
an Eulerian grid is applied by iLivnd (11993 ). Most explicit Eulerian approa ches 



were based on versions of the Prometheus implementation (Fr vxell et al.L 11989 ) of 



I ap 



the Piecewise ParaboHc Method (PPM) proposed b y jColella & Woodward ( 1984t) 
These include the simulations of Arnett & Livnel (i l994ab). Khokhlov (1995) 
Niem ever & Hillebrandd ( Il995bh: INiemever et alj (11996.). .Reinecke et al., 



Hillebrandt et al 



(I2000I). 



Reinecke et al. 

r-r 



;r et al J (11996.). .Reinecke et al. (■1999a) 
(l2002bh. iGamezoetalJ (|2003[J2004|) 



Ropke & Hi llebrandt (2004), Calder et al. (2004), Plewa et al. (2004), Ropke (2005) 

Ropke & Hillebrandt (.2005a.b.) . .Ropke et al.. (.2005a. .2006a) an d .Gamezo et al.l 

(l2005h . An alternative approach was taken by Bell et al.l (l2004ah . who modified 



This 



the hydrodynamical equations to account for low-Mach number flows only, 
provides an efficient scheme to numerically simulate subsonic deflagrations. 

The facts that the WD star is expanding and that the flame is not necessary 
domain-filling in all stages of the explo sion is taken in t o account in th e grid-based 
schemes in various approaches. While iReinecke et al.l (Il999al l2002bl) use a static 
computational grid with a fine-resolved central part and a coarse outer grid to ac- 
count for the expansion, iRopkd (l2005b implemented a uniform moving grid that 
co-expands with the exploding WD. Since on average the explosion process is more 
or less spherical, one can gain resolution in the part of the star occupied by the flame 
with nested moving grids ( IRopke et al.l l2006a). A different technique providing res- 
olution where needed is adaptive mesh refin ement as applied in the simulations by 
Gamezo et all (12003) and lCalder et"aD (l2004 . 



3.4.2 Flame model 



In earlier works (IMiiller & Arnetti Il986l: iLivnel Il993h the flame was advanced ac- 
cording to discrete boolean criteria. A computational cell was burned if certain 
conditions in neighboring cells were fulfilled. This strategy is clearly dependent 
on the grid geometry and therefore in recent simulations different approaches have 
been taken. Two major strategies to tackle the pro blem of the unresolved inter- 



nal flame structure can be distinguished. Khokhlov (i 19931 11994 11995) adapted a 



flame capturing technique that mimics flame propagation by an artificial diffusion 
mechanism which broadens the internal flame structure to a certain number of com- 
pu tational grid cells. This me t hod was applied in t he SN la e xplosion simulation s 
by iGamezo et all (l2003l Eooi l2005l) . ICalder et alJ (l2004. andlpTewa et all (l2004 . 
A completely different approach was taken by Reinecke et al.l (Il999bl) . who treat 
the flame as a sharp discontinuity separating the f uel from the ashes. It is numeri- 
cally represented applying the level-set technique (lOsher & Sethianlll988tlSethian , 
19961) . Here the flame front is associated with the zero level set of a scalar function 
G representing the distance from the interface. A model f or flam e propagation based 
on this technique was developed by ISmiljanovski et al.l (Il997h and a modified ver- 
sion for thermonuclear flames in SN la explosion was presented by iReinecke et al 
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( Il999t » (for details of the implementation see also Hillebrandt et al., 2005!). This 



scheme was applied in the simulat i ons by Hillebrandt et al. (2000), Reinecke et al. 



([2002c a.b ). iRopke & HiUebrandtl (1200 4. 2005b a), Ropke (2005), I Schmidtetal. 



j2005q) , iSchmidt & Niemeveil (l2006l) , and iRopke et al.i(i2005aii2006: 



3.4.3 Turbulent combustion model 

As outlined above, SN la explosions are believed to at least start out subsonically 
in the deflagration mode of flame propagation. Thus the flame is subject to the 
interaction with turbulence generated by instabilities. This mechanism turns out 
to be crucial to accelerate the flame propagation and must therefore be taken into 
account in any valid SN la model. 

The theory of turbulent combustion has been extensively developed for terrestrial 
combustion phenomena (see Peters. 2000) . The nuclear reaction kinetics in ther- 
monuclear combustion waves is much simpler than the chemical reactions in terres- 
trial flames and due to the high thermal conductivity of the degenerate electron gas in 
WD matter some of the characteristic nondimensional numbers differ Nonetheless, 
the turbulent combustion process in the deflagration mode in SNe la bears striking 
similarity to premixed turbulent flames (i.e. turbulent flames in mixtures of fuel and 
oxidizer) in terrestrial combustion processes. 

The wide range of scales involved in turbulent combustion phenomena renders 
direct simulations impossible for most situations. Therefore only parts of the inter- 
action range of the flame with turbulence and the resulting surface enlargement of 
the flame can be resolved. This is usually compensated by attributing an effective 
turbulent flame speed Sx to the unresolved flame front, which must be determined by 
theoretical considerations. One of the cornerstones of the theoretical description of 
turb ulent combustion is the notio n of different regimes of flame/turbulence interac- 
tion (Niem ever & Kersteinlll997h . These regimes are distinguished by the ability of 



turbulent eddies to penetrate the internal flame structure. Since the Gibson scale is 
much larger than the flame width for most parts of the SN la explosion, this will not 
be the case here and accordingly the combustion falls into the regime of wrinkled and 
corrugated flamelets. Here, the full flame structure is corrugated by the interaction 
with turbule nce and the resulting surface enlargement accelerates its propagation. As 



first noted by Damk ohleri (Il940h . the flame propagation in this regime completely de- 



couples from the microphysics of the burning for sufficiently strong turbulence. It is 
entirely determined by the turbulent velocity fluctuations, that is, it is proportional 
to the turbulent velocity fluctuations. 

One of the challenges of deflagration models of SN la explosions is thus 
to determine these velocity fluctuations correctly. Since the resolution in multi- 
dimensional simulations is insufficient to resolve the phenomena directly, modeling 
approaches have to be taken. In the simplest models, a minimal propagation ve- 
locity o f the flame is artificial l y imp osed, which is not p hysically determined in the 
model ( Ga rcia-Senz & Bravol l2005h . Contrary to this, lOamezo et al.l (|2003), and 



Calder et al. (2004) assume the flame to be driven by buoyancy-induced instabilities 



on unresolved scales giving rise to a turbulent flame speed of Si - 0.5 ^jAgl, where 
g and / denote the gravitational acceleration and the computational grid cell size, re- 
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spectively. A is the Atwood number derived from the density contrast over the flame 
front. A more sophisticated approach was proposed by Niemeyer & Hillebrandt 



( Il995ah . Guided by the techn ique of Lar^e E ddy Simulations they implemented a 
turbulent subgrid-scale model ( IClementlll993h in SN la explosion simulations. This 
model determines the turbulence energy on unresolved scales based on conservation 
laws. The hyperbolic nature of the equations of hydrodynamics, however, does not 
allow for a closed system of equations here so that closure assumptions have to be 
invoked. A physical ly better motivated ap proach based on localized closures was 
recently proposed bv lSchmidt et alj (l2005bllcl) . 

In the very late stages of the SN la explosion the fuel density drops due to expan- 
sion of the WD to values where the flame width becomes broader than the Gibson 
length. Then, turbulence penetrat es the internal structure of th e flame and it enters 
the regime of distributed burning f Niemever & Kerstein, 1997). Here, diff^erent scal- 
ing laws for the turbulent flame speed apply ( Damkohler, 1940) and a prelim inary 
test of the effects on SN la models was given bv lRopke & Hillebrandtl (l2005bl) . 



3.4.4 Nuclear reactions 

The computational expenses for a full nuclear reaction network are prohibitive to 
run it concurrently with the explosion hydrodynamics simulations in present multi- 
dimensional models. Since here the only dynamically relevant parameter is the en- 
ergy release in the reactions, it is sufficient to apply a simpliffed description. Usually, 
only a few nuclei representative for the fuel mixture, the iron group elements and the 
intermediate mass elements in the ashes and effective reactions between them are 
accounted for (e.g. Reinecke et al.L '2002aV 

However, in order to derive observables (such as lightcurves and spectra) from 
the explosion models, the exact chemical composition of the ejecta needs to be 
known. This can be achieved by advecting a number of tracer particles in the ex- 
plosion hydrodynamics simulations which record the evolution of temperature and 
density. This adds a Lagrangean component to the Eulerian code. Based on the 
data gained from the tracer particles it is then possible to a posteriori reconstruct the 
detailed nuclear reactions in the burnt material with extended nuclear reaction net- 



works. This technique was implemented for SN e la by iTravaglio et alJ (120041) and 
also applied in the study bv lRopke et al.l (12005 al) . 



3.4.5 Results of multi-dimensional models 

Whi le earlier two-dim ensional (2D) deffagration models of SN la f ailed to explode 



(e.g. [ Khokhlovlll995h others succeeded to (weakly) unbind the WD lNiemeyer et al. 
(119961) . Nowadays, it seems that a consensus has been reached on the gen 



eral capability of the 30 version of this mo d el to give rise to robust expl osions 
dReinecke et all l2002bt iGamezo et all 120031; iRopke & HiUebrandtl l2005al) . Al- 
though the energetics (asymptotic kinetic energy up to 7 x 10^° erg) and the pro- 
duction of ^''Ni (~0.4Mq) as the main global characteristcs fall into the range of 
observational expectations (iContardo et al.L uOOOt IStritzinger et al.L l2005h . they are 
still on the weaker side of "normal" SNe la. 
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Since the initial flame configuration in the explosion is undetermined yet (see 
Sect. 12. II ). several studies addressed its eff'ect on the outcome of the simulations. Off"- 
center ignitions were analyzed based on 2D - simula tions by iNiemever et al.l (119961) 
and 3D-full star simulations by ICalder et al.l (12004 V The latter study indicated that 
a single perfectly spherical ofi'-center initial flame quickly em erges to the s urface 
of the WD and fails to burn sufficient material to explode it. iPlewa et al.l (120041) 
suggested that the material breaking through the surface of the still gravitation- 
ally bound WD may c ollimate on the opposite sid e of it, evoking a detonation 
there. Contraty to that, iRopke & Hillebrandtl (l2005al) found that a more structured 
off^-center initial flame (motivated by pre-ignition convection, cf. Sect. 12. Il l stays 
more or less in place and gives rise to a viable explosion. The effects of the num- 



ber of ignition spots have been addressed in s everal 2D studies (INiemever et al. , 
1996t iReinecke et al.i Il999al: iLivne et al.L l2005b giving a wide range of results de- 
pendent on the number a nd distr ibution of the initial flame kernels and in 3D sim- 
ulations jReinecke et all 12 002b': Ropke & Hillebrandtl l2005al: iRopke et all l2006at 
Garcia- Senz & Bra vol |2005|). These simulations indicate that multi-spot ignition 
scenarios can lead to more vigorous explosions than centrally ignited setups and 
put the deflagration scenario in better agreement with observations. Interestingly, 
the 3D simulations, although also exhibiting a depencence on the number of igni- 
tion spots, gave much more robust results than the 2D studies (Ropke et al., 2006a). 
Similar conclusions were drawn on the basis of a stochastic ignition model in 
which a spatial and temporal distribution of the flame ignition spots was assumed 
ISchmidt & Niemeyer, 2006). 

A successful explosion, however, is not the only requirement for a valid SN la 
model. Observables derived from such models must match the observations of 
SNe la. Moreover, the model must allow for a certain variation of the results in 
order to reproduce the observed range of diversity in the characteristics. No general 
consensus has been reached so far regarding these questions. Besides the low ex- 
plosion energies and ^^Ni production the problems of most deflagration simulations 
include unburnt material left behind at low velocities in the ejecta due to downdrafts 
in the buoyancy-induced large-scale flame pattern and low production of intermedi- 
ate mass elements. Recently, a series of 3D simulaions was analyzed to determine 
wheter these shortcomings are generic to the deflagration model or caused by the 
simplicity of some setups. It seems that the model is generally capable of reproduc- 
ing main observational features, but some of the problems persist (see Sect. |4] for 
details). 



Gamezo et al.l ( 120041) claim that a delayed detonation is necessary to achieve 
agreement with observations. In their 3D model, a detonation front is artificially 
initiated after a deflagration stage at a pre-selected time and location. As ex- 
pected, it travels through the WD burning most of the materinal in the inner parts. 
A 2D implementation of the d elayed detonation model was recently presented by 
Golombek & Niemeven (120051) . Here, the DDT is assumed to occur when the defla- 
gration flame enters the distributed burning regime which appatently happens first in 
the outer parts of the flame front. This rises the question of whether the detonation 
(if not allowed to travel through ashes) can catch up with the expanding material 
on the opposite side of the WD. 2D simulations by ISteinmetz et al.l (119921) studied 
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pure prompt detonation in rapidly rotating WDs. Although here regions may exist 
with densities allowing for the production of intermediate mass elements, the authors 
conclude, that the ratio of iron group elements to intermediate mass elements arising 
from such models is inconsistent with normal SNe la. 

To settle the question of the explosion mechanism, detailed analysis of the dif- 
fernt models and comparison with observations on the basis of synthetic light curves 
and spectra are required. 



3.5 Complementary small-scale simulations 

Complementary studies of flame propagation on small scales focus on effects in only 
a narrow window in scale space. These therefore reach a much higher resolution and 
can model the processes in a more realistic way than SN la simulations on scales of 
the WD star Consequently, they give insight into mechanisms that are unresolved 
there either testing the validity of assumptions or providing input data to these mod- 
els. 



On the basis of resolved ID flame simulations. lTimmes & Wooslevi (119921) deter- 
mined the laminar speed of deflagration flames in SNe la and lSharpd (11999) analyzed 
the propaga tion velocity of detonations. Both serve as input in m ulti-dimensional 
simulations (iReinecke et al.[l2002btlGolombek & Niemeyenl2005h . 

Other studies concerned the validity of subgrid-scale models of the interaction 
of the fl ame with turbulence (Schmidt et al.Ll2005 a). This interaction was also simu- 
lated bv lNiemever et alj (Il999h . The deflagration flame propagation on small scales 



in SNe la subject to the Rayleigh-Taylor i nstability was tested in 2D-simulations 



m tiJNe la subject to tne Kayieign- l ayior 

by Ehokhiovl / l993l) an d iBeU et alJ (l2004ch and in 3D-setups by iKhokhlovl (119951) 
and Zingale et alj (l2005b . The effects of the Landau-Darrie us instabihty on the ther- 
monucl e ar flame have been analyzed i n 2D-simulations by iNiemever & H illebrandtl 
(Il995al). iBlinnikov & Sasorovl d 19961). iRopke et al.l (l2003h . iRopke et all (i2004a.b,') . 



and Bell et al 



(l2004bt) . 



Dursi et al.l (l2003h studied the response of thermonuclear 



flame to curvature and stretch. Most of these simulations confirm the assumptions 
made in large-scale SN la models. 

The question whether a transition of the turbulent burning regime from flamelet 
to distribut ed burning would cause a deflagration-to-detonation transition was ad- 
dressed by iLisewski et al.l (l2000blah . who conclude that such an event is very un- 
likely. 



4 Example: A deflagration SN la model 

To illustrate the numerical modeling of SN la explosions in three spatial dimensions, 
a deflagration model in a specific implementation shall be discussed in this section. 
It starts out with a cold isothermal (T - 5x 10^ K) Chandrasekhar-mass WD. Unless 
otherwise stated, its composition was assumed to consist of equal parts of carbon 
and oxygen. 
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4.1 Numerical techniques 



The numerical techniques forming the foundation of the im plementation of the defla- 
gration SN la model presen ted here are described by Reinecke et al.l( 1999br) . lRopke 
(l2005h . ISchmidt et all (|2005t|), and |Niemeiier& Hillebrandt (1995^). The funda- 



mental approach is that of Large Eddy Simulations, where hydrodynamics on the 
resolved scales is modeled in ?i fini te volume approach based on the Prometheus 
implementation dFrvxell et al.Lll989l) . 

The equation of state describing WD material contains contributions from an 
arbitrarily degerneate and relativistic electron gas, an ideal gas of nuclei, radiation, 
and possibly electron-positron pair creation and annihilation. 

As a consequence of|2]of the enumeration in Sect. 13. II a subgrid-scale model is 
applied to account for turbulent effects o n unresolved scales. Som e older simulations 
follow the implementation suggested by lNiemeyer & Hillebrandtl(ll99 5b). wh ile one 



recent highly resolved run applied the updated modeling approach of ISchmidt et al. 
(l2Q05b.) . 

According to [T] of the enumeration, the flame representation is achieved in a 
modeling approach. Seen from the scales of the WD star, the flame appears as a 
sharp discontinuity separating the fuel from the ashes. A suitable numerical method 
to follow the evolution of such an in t erface i s the level set technique introduced 
to SN la modeling by iReinecke et al.l (Il999al) . In this implementation, the flame 
propagation speed needs to be prescribed. This quantity, however, is not arbitrary in 
the presented modeling framework. For turbulent combustion in the flamelet regime 
(cf. Sect. I3.4.3I I, which applies to the burning in major parts of SN la explosions, it 
is given by the turbulent velocity fluctuations. These are determined by the subgrid- 
scale turbulence model. 

Nuclear reactions are implemen ted in the simplified approach (cf. Sect. 13.4.41 ) 
proposed by IReinecke et al.l ( l2002ah . The progenitor material is composed of '^C 
and '^O. At high fuel densities nuclear burning terminates in nuclear statistical equi- 
librium represented by a mixture of ^*Ni and a-particles. Once the fuel density drops 
below 5.25 X 10^ g cm"-', burning will become incomplete and intermediate mass el- 
ements (represented by ^^Mg) are produced. The respective difference in nuclear 
binding energy is released which provides sufficient accuracy to model the dynamics 
of the ex plosion. The chemical composition of the ejecta is derived in a postprocess- 
ing step (ITravaglio et al.ll2004l) . 



4.2 Results of Simulations 

Several simulations based on the implementation described abov e , both in two a nd 
three spatial dimensions, have been presented by Reinec ke et al.l (Il999a[ l2002allbl) . 
In the 2D simulations, numerical convergence in the global quantities was demon- 
strated. For the i mplementation on a co-expanding computational grid, a similar 
result was found by Ropkg (l2005h (see Fig. [T]!. The numerical convergence naturally 
arises from the interplay of the resolved flame front representation with the turbulent 
subgrid-scale model. Ideally, a lack of resolution of large-scale features in the flame 
front representation should be compensated by an increased turbulent flame prop- 
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Figure 1 : Total energy in 2D-simulations with co-expanding grid for different num- 
bers of computational grid cells. 



agation velocity determined from the subgrid-scale approach. Of course, a certain 
threshold of resolution will need to be exceeded to reach this regime in the numerical 
implementation. 

One requirement to reliably derive observables from simulations is that the evo- 
lution of the models be followed to the stage of homologous expansion. In this 
hydrodynamically relaxed situation, the velocity of the ejected material is propor- 
tional to its radius. Obviously, a static computational grid that sufficiently resolves 
the flame propagation does not allow to follow the expansion over sufficiently long 
time scales. Therefore iRopkd (1200 5^ proposed to use a moving computational grid 
that co-expands with the WD. With this implementation it is in principle possible to 
follow the evolution for arbitrary times, but JRopkel 120051) showed that simulating 
the first 10 s after ignition is sufficient to reach homology with reasonable accuracy. 
An example of such a simulation is shown in Fig.|2] where the isosurface represents 
the zero level set of G and is associated with the flame front. This simulation was 
carried out on only one octant of the WD assuming mirror symmetry to the other oc- 
tants. The flame was ignited centrally with a toroidal perturbation superposed to the 
spherical shape. Of course, at times after ~2 s when burning terminates in the model, 
the zero level set of G looses physical meaning, but it still indicates the approximate 
boundary between unburnt material and ashes. 

Due to the high computational expenses, most 3D simulations (such as the one 
described above) comprise only one octant of the WD. However, only full-star setups 
allow to account for asymmetr y effects. On the basis of spectrapolarimetry obser- 
vations of several SN e la (e.g. Wang et all 2003h t hese are expected to occur in at 
least some explosions. iRopke & Hillebrandtl (12005 ah showed, that in the deflagration 
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Figure 2: Snapshots from a single-octant SN la simulation performed on a uni- 
form expanding computational grid. The isosurface corresponds to G = and is 
color-coded with a measure of the turbulence strength derived from the subgrid-scale 
model (top left to bottom right: ignition flame, 0.5 s, 2.0 s, and 10.0 s after ignition). 



model such asymmetries arise exclusively from irregularities in the flame ignition 
conditions and not from large-scale instabilities and resulting prefered modes in the 
flow patterns. Therefore, fixing a symmetric initial flame shape and studying the 
influence of other physical parameters on single-octant explosion models is a valid 
approach (see Sect. l4.4l below). 

To illustrate the typical flame evolution in defl agration SN la models, the full-star 
model presented bvl Ropke & Hillebrandti (12005 al) shall be described here. It started 
out from an asymmetric initial flame configuration which was set up by randomly 
distributing spherical flame kernels around the center of the WD. This resulted in a 
foamy structure slightly misaligned with the center of the WD. 

Starting from this initial flame configuration (shown in Fig. [3]), the evolution of 
the flame front in the explosion process is illustrated by snapshots of the G = iso- 
surface at f = 0.3 s and t = 0.6 s in Fig. [3] The development of the flame shape from 
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Figure 3: Snapshots from a full-star SN la simulation starting from a multi-spot 
ignition scenario. The logarithm of the density is volume rendered indicating the 
extend of the WD star and the isosurface corresponds to the thermonuclear flame. 
The last snapshot corresponds to the end of the simulation and is not on scale with 
the earlier snapshots. 



ignition to f = 0.3 s is characterized by the formation of the well-known "mushroom- 
like" structures resulting from buoyancy. This is especially well visible for the bub- 
bles that were detached from the bulk of the initial flame. But also the perturbed parts 
of the contiguous flame closer to the center develop nonlinear Rayleigh-Taylor like 
features. During the following flame evolution, inner structures of smaller scales 
catch up with the outer "mushrooms" and the initially separated structures merge 
forming a more closed configuration (see snapshot at f = 0.6 s of Fig. O. This is 
a result of the large-scale flame advection in the turbulent flow and the expansion 
of the ashes. Up to this stage the flame was strongly anisotropic. However, in the 
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later evolution a preferentially lateral growth of bubbles filled with ashes smoothes 
out parts of the anisotropies. The flame develops a more spherical shape and only 
a slight anisotropy is retained. After about 2 s self -propagation of the flame due to 
burning has terminated in the model. The subsequent evolution is characterized by 
the approach to homologous expansion. The resulting density structure at the end of 
the simulation is shown in the f = 10 s snapshot of Fig. [3] The flame evolution agrees 
with the expectations outlined in Sect. 12.21 

Since the flame ignition process is not modeled in the explosion simulations, 
there exists considerable freedom in choosing the number and distribution of igni- 
tion kernels (cf. Sect. 12.11 ). The results from models starting with a central ignition 
(cf. Fig. |2| indicate that such a setup leads only to weak explosions. Recently, the 
capabilities of multi-spot ignition models have been analyzed in a systematic way. 
This was possible due to a modification of the moving grid implementation to two 
nested grids. An inner, fine-resolved grid f ollows the flame pro pagation while an 
outer coarse grid tracks the WD expansion JRopke et al.L l2006ah . In this way, it is 



possible to accumulate a large fraction of the available computational cells in the 
inner regions thereby resolving detailed multi-ignition configurations. Confining the 
ignition volume around the center of the WD, it may be expected that there exists a 
certain number of ignition spots which maximizes the burning and energy release. 
A sparse ignition would decrease the flame surface and thus the burning rate, while 
a too dense distribution of flame kernels will lead to a rapid merging of flame parts 
again decreasing the surface area. This was confirmed by the simulations, which 
however did not sharply single out an optimal number of ignition points but revealed 
a rather robust behavior over a wide range. Extreme cases, however, reproduced the 
anticipated dependence. 

4.3 Comparison with observations 

Due to recent progress in deriving observables from multi-dimensional deflagration 
simulations, a direct comparison with details of observations of nearby SNe la has 
come into reach. Since these contain no other parameters than the initial conditions, 
the question arises if the outcome of simulations as the ones described above meets 
the observational constraints. Such constraints result from the global characteristics 
derived from observations, observed light curves, and spectra taken from nearby SNe 
la. 

The global characteristics derived from SN la observations state that a valid ex- 
plosion model should release arou nd 10^^ erg of en e rgy and produce ~0.4 . . . 0.7 M p 



of ^^Ni in the nuclear burning dContardo et al.l uOOOt IStritzinger et al.l l2005b . 



However, there exists a large diversity in the observations ranging from the class 
of sub-luminous SNe la (like SN 1991bg with probably -0.1 Mq of ^^Ni) to 
super-luminous events (e.g. SN 1991T with a ""^Ni mass close to 1 Mq). Deflagration 
models started with a multi-spot ignition setup typically possess 6 ... 8 x 10^*' erg of 
asymptotic kinetic energy of the ejecta. They produce up to ~ 0.4 Mq of ^''Ni. Thus 
they fall into the range of observational expectations, but in the current stage do 
not account for the more energetic SNe la. One reason for the low energetics may 
be that in the presented implementation nuclear burning is assumed to cease when 
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Figure 4: Synthetic light curves derived from model 2_2_2 of (IRopke et all 12005a) 
(solid curves) compared with observed light curves from SN 1994D. 



the fuel density drops below 10^ g cm ^, because the flame is then expected to enter 
the distributed burning regim e. An approach to handling this stage was given by 
Ropke & HillebrandtllEoOSbl) . 



4.3.1 Lightcurves 

A further requirement is that synthetic lightcurves agree with observed ones. These 
are sensitive to the energy release, the ^^Ni production, as well as to the distribution 
of elements in the e jecta . In Fig . |4] synthetic light curves derived from the 2_2_2 
simulation (Ropke et al.1 l2005a|) are compared with observ ations of SN 1994D' 



efail 
)05rst 



(IBlinnikov et alil2005l see also lSorokina & Blinnikovi 12003). The multi-band light 



Blinnikovet al.1 (11993) 



cur ve of this model was calculate d using the Stella code of 
and IBlinnikov & Sorokiiial (|2000|) . 

The model produced 0.3 Mq of ^^Ni, although observations require somewhat 
higher ^^Ni mass (~ 0.4 Mq) for the assumed distance. Nonetheless, there is gener- 
ally very good agreement in the B and V bands near peak luminosities and in decline 
rate 20 days after the peak which is most important for cosmological appUcations of 
type la supernovae. 



*The distance to SN 19 94D is still controversial and liere tlie value of 30.4 is adopted from 
iDrenkhahn & Richtlej Il999h 
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Figure 5: Synth etic nebular spectrum compared with observations (from 
Kozmaetalll2005h . 



4.3.2 Spectra 

A much harder test for the models is posed by the comparison of synthetic with 
observed spectra since these depend on details in the composition of the ejected 
material. 

Nebular spectra provide a means of studying the central parts of the ejecta, since 
they are taken at epochs where these have become transparent due to expansion. 
Thus they explore the "heart of the explosion" and are a valuable tool to study de- 
tails of the physical processes involved in the explosion stage. Unfortunately, only 
one single synthetic late time spectrum is available from deflagration SN la models 
( IKozma et al.l uOOSb . It was derived from the simulation shown in Fig.|2] The arti- 
ficial and simple initial flame shape chosen here gives reason to not expect a good 
agreement between model and observation. However, as shown in Fig.|5] the broad 
iron features of the observed spectra are qualitatively reproduced. An inconsistency 
of the model with the observed nebular spectra is the appearance of a pronounced 
oxygen line at 6300 A. Both features of the synthetic spectrum share a common 
origin. The broad iron lines are caused by NSE material that is transported in the 
uprising plumes of ashes and thus distributed in velocity space. At the same time 
strong downdrafts carry unburnt material towards the center of the WD. 

The disagreement may in part be attributed to the simplicity of the explosion 
model. Its highly symmetric initial flame shape with large imprinted perturbations 
favors a pronounced evolution of large-scale Rayleigh-Taylor features. This prob- 
lem, however, is not necessarily generic to all SN la deflagration models. Multi-spot 
ignitions may lead to a more complete burning in the central parts of the WD. This 
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Table 1 : Variation of initial parameters in SN la explosion models. 



Parameter 


range of variation 


effect on ^^Ni 


effect on total 






production 


energy 


X(12C) 


[0.30,0.62] 


<2% 


-14% 


Pe[109g/cm3] 


[1.0,2.6] 


-6% 


-17% 


Z[Zo] 


[0.5,3.0] 


-20% 


none 



was demonstrated bv lRopke etal.l ( l2006ah . For an optimal number of ignition bub- 
bles (and a rather wide range around this number), the central parts of the ejecta 
become dominated by iron group elements. Whether the achievable suppression 
of unburnt material at low velocities is sufficient to be consistent with observations 
needs further study. 

A powerful diagnostic tool to compare SN la models with observations is pro- 
vided by the abundance tomography presented by lStehle et al.l (120051) . It makes use 
of spectra taken from SN 2002bo with an extraordinary good time coverage. Fitting 
this sequence of data with synthetic spectra unveils the composition of the ejecta in 
velocity space slice by slice, since the photosphere moves gradually inwards with 
the expansion of the remnant. This abundance tomography of the ejecta can be com- 
pared with results of 3D models, when av eraged over the ang les. Qualitatively, the 
mixed composition of the ejecta found by lStehle et al.l (l2005h is reproduced by de- 
flagration SN la models in a natural way since these predict a distribution of burnt 
material with the rising bubbles. A problem was, however, that older predicted 
large unburnt mat erial fractions in th e central parts of the ejecta in disagreement 
with the results of S tehle et al.l (1200 5). A recent high-resolved simulation cures this 
problem by clearly re producing the iron-group dominance in the low-velocity ejecta 
(lRopkeetalll2006bl) . 



4.4 Diversity and correlations 

The recent developments in the deflagration SN la explosion modeling outlined in 
the previous section seem to indicate that such a model is capable of reproducing the 
main features of observed objects but they do not rule out the alternative of a delayed 
detonation. 

Here, the question of how multi-dimensional SN la models can be applied to 
tackle questions from the cosmological applications of these objects shall be illus- 
trated on a set of simple deflagration models. Of particular interest is how such mod- 
els behave under variation of physical parameters. Do they reproduce the observed 
diversity of SNe la? Are correlations between observables evident in the results? 

Unfortunately, 3D deflagration models of SNe la as described above are com- 
putationally expensive. To moderate these expenses, simplified setups may be used 
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Figure 6: Snapshots of the flame front evolution at f = 1.0 s after ignition from 
models with p^ = 1.0 x 10'' g cm^^ and ^('^C) = 0.42 (left); pc = 2.6 x lO** g cm^^ 
andX(i2c) = 0.42 (middle); p,. = 2.6 x lO^gcm^^ wdXC^C) = 0.62 (right) 



to study eff'ects of ph ysical parameters on the explosion models. Such an approach 
was recently taken by lRopke et al.l (12005 al) and resulted in the first systematic study 
of progenitor parameters in 3D models. The basis of this study was a single-octant 
setup with moderate (yet numerically converged) resolution. However, the lack in 
resolution did not allow a reasonable multi-spot ignition scenario and thus only weak 
explosions can be expected. It was therefore not possible to set the absolute scale of 
effects in this approach, but trends can clearly be identified. 

The parameters chosen for the study were the WD's carbon-to-oxygen ratio, its 
central density at ignition and its ^^Ne mass fraction resulting from the metallicity 
of the progenitor All parameters were varied independently to study the individual 
effects on the explosion process. In a realistic scenario, however, these parameters 
are interrelated by the evolution of the progenitor binary system. The results of this 
survey are given in Tab. [T] 

A variation of the carbon-to-oxygen ratio affects the energy production in the 
burning due to the differences in the binding energies of these two nuclei. Counter- 
intuitively, this results in no significant change in the ^*Ni production. This finding 
is consistent with the flame very similar flam e evolution in simulations wi th different 
carbon mass fraction of the WD (cf. Fig. |6]l. iRopke & Hillebrandti (120041) explained 
this effect by the fact that the potentially higher energy release in carbon-rich mod- 
els is buffered by a higher a-particle fraction in the NSE material and only released 
when burning is already incomplete. 

Models with lower central densities show a delayed flame evolution (see Fig.|6ll 
and consequently a lower and delayed energy production. This is due to the fact that 
the flame experiences a lower gravitational acceleration in these models resulting in 
decreased turbulence generation. Therefore less ^^Ni is produced in these models. 
This effect is even more pronounced due to the fact that in low-density WDs less ma- 
terial is present at sufficient densities to be potentially burnt to NSE. A counteractive 
effect is expected at higher densities. Here, electron captures will become important 
favoring neutron-rich isotopes over ^''Ni in the NSE. The dynamical effects of elec- 
tron captures are, however , not yet implemented in the explosion model and therefore 
the survey of iRopke et al.l (l2005al) is restricted to relatively low central densities. 
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Figure 7: Peak luminosity vs. decline rate of the light curve in the B band (di- 
a monds correspo nd to SN la explosion models). Compared with original relation 
bv lPhillips et alj (11999.) (dashed curve) and shifted relation (solid curve) in the left 
panel. 



A higher metallicity of the main sequence progenitor star results in an increased 
^^Ne mass fraction in the WD. This is an isotope with neutron excess and therefore 
ag ain favors the produc tion of neutron-rich species over ^'' Ni in the NSE. The re sults 
of iRopke et alj (l2005al) confirm t he analytic prediction bv lTimmesetal.1 (l2003h and 
agree with lTravaglio et al.l (l2005h . The metallicity parameter, however, has no effect 
on the explosion dynamics and the energy production in the models. 

To determine the efl'ects of these variations on observables, synthetic light curves 
were derived from all models (an example is shown in Fig. |4]l. From these, the peak 
luminosities and decline rates (in m agnitude s 15 da ys after maximum; Amis), were 
determined. The pioneering work bv lPhillipsI d 1993b estabUshed the relation between 
the peak luminosity and Amis as one of the primary tools to calibrate cosmological 
SN la distance measurements. The so-called Phillips relation quantifies the decrease 
of Ami5 for brighter SNe la. 

The results from the presented deflagration models are compared with the rela- 
tio n given by Phillips et al. (1999) in Fig. [T] Obviously, the absolute magnitude of 
the lPhillips et al.l (Il999h relation is not met by the set of models (cf. the upper panel 
of Fig. |7]). Moreover, the range of scatter in Amis is much narrower than that of 
the set of observations used by Phillips et al. (1999). The simulation s exhibit a large 
scatter but are consistent with the slope of the IPhilhps et al.l (119991) relation. It is, 
however, obvious, that a better agreement cannot be expected given the fact, that the 
parameters in the survey were chosen independently. A consistent stellar evolution 
would pick a sub-sample of the set of models possibly narrowing the range of scatter 

Nonetheless, with the set of models, the question can be answered, which param- 
eter dominates the slope in the direction of dimmer events for faster decline rates. 
The varied parameters are coded by different line styles in the right panel of Fig. |7] 
Clearly, the progenitor's metallicity can be identified as this parameter. Variations in 
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the central density and the carbon mass fraction of the WD superpose a scatter on 
the dominant relation. 



5 Conclusions 

The rapid development of multi-dimensional SN la models over the past years led 
to an better understanding of the mechanism of these events. In particular, it has 
been shown, that 3D SN la simulations on the basis of the deflagration model are 
capable of reproducing main features of observed SNe la. However, currently it 
seems that they can only account for the weaker objects. A better determination of 
the initial conditions, i.e. the progenitor evolution and flame ignition, may improve 
the situation. It may, however, also turn out that the current model is incomplete. 
A delayed detonation was suggested by several authors to improve the agreement 
with observations, lacking however a convincing mechanism for a deflagration-to- 
detonation. 

On the exemplary case of a particular numerical implementation of the deflagra- 
tion scenario it was shown how such a model can be validated against observations 
and applied to determine the origin of the diversity of SNe la as a first step to theo- 
retically assess the empirical relations utilized in SN la cosmology. 
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